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ABSTRACT 



We use Swift/BAT Earth occultation data at different geomagnetic latitudes 
to derive a sensitive measurement of tlie Cosmic X-ray background (CXB) and 
of tlie Eartli albedo emission in the 15-200 keV band. We compare our CXB 
spectrum with recent [INTEGRAL, BeppoSAX) and past results (HEAO-1) and 
find good agreement. Using an independent measurement of the CXB spectrum 
we are able to confirm our results. This study shows that the BAT CXB spectrum 
has a normalization ~ 8 ± 3 % larger than the HEAO-1 measurement. The 
BAT accurate Earth albedo spectrum can be used to predict the level of photon 
background for satellites in low Earth and mid inclination orbits. 



Subject headings: cosmology: observations - diffuse radiation - galaxies: active 
X-rays: diffuse background - galaxies - Earth 



1. Introduction 



There is a gen eral consensus that t he cosmic X-ray background (CXB), discovered more 
than 40 years ago (iGiacconi et al.lll962l ). is produced by integrated emission of extra-galactic 



point sources. The deepest X-ray surveys to date (IGiacconi et al.l l2002l : [Alexander et al. 



20031 : iHasingerl 120041 ) have shown that up to virtually 100% of the < 2 keV CXB radi- 



ation is accounted for by Active Galactic Nuclei (AGN) hosting accreting super-massive 
black holes (SMBHs). However, the frac tion of CXB emissio n resolved into AGNs declines 
with energy being < 50% above 6keV (IWorsley et al.l 120051 ). The unresolved component 
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may be attributed to the emission of a yet undetected population of highly-absorbed AGN. 
These AGN should be characterized by havin g column densities ~ 10^"^ H-atoms cm~^ and 
a space density peaking at redshift below 1 (IWorsley et al.ll2005l ). Such a population of 



Compton-thick AGN is invoked by po pulation synthesis models (e.g. IComastri et al.l Il995 



reister fc Urryl 2005 : iGilli et al.ll2007l ) to reproduce the peak of the CXB emission at 30 keV 



flMarshall et al.lll980l ). 



Thus, an accurate measurement of the CXB spectrum in the 15-200 keV energy range 
is important to assess and constrain the number density of Compton-thick AGNs. Such 
measurements are complicated by the fact that instruments sensitive in this energy range 
are dominated by internal detector background and are not designed to measure the CXB 
spectrum directly (excluding HEAO-1 A2). The typical approach is to produce an ON — 
OFF measurement, where taking the difference between the ON and the OFF pointings 
eliminates the internal background component. 

There are different methods to obtain a suitable OFF observation; the HEAO-1 mea- 
surement of the CXB spectrum in the 13-180 keV range was obtained by blocking the aper- 
ture with a movable Csl crystal. Also the Earth disk can be used to modulate the CXB 
emission. This approach is the o ne used in recent CXB intensity measurem ents performed 
by INTEGRAL and BeppoSAX Jchurazov et ahlboO?! : IProntera et aPboO?! ). 



Here we report on two independent measurements of the CXB emission using Swift/BAT. 
For the first method, we use the Earth occultation technique similarly to the INTEGRAL 
and BeppoSAX analyses while for the second one we make use of the spatial distribution of 
the BAT background. 

The structure of the paper is as follows. In §[2] we present the details of the observations 
and describe the BAT background components. We also derive a rate-rigidity relation which 
is fundamental for suppressing the background variability due to Cosmic Rays (CRs). In §[3] 
we present the details of the Earth's occultation episodes undergone by BAT and the analysis 
method for the occultation measurement. In § HI we discuss all sources of uncertainties 
which affect our occultation measurement which is then presented in § O The alternative 
measurement used to verify the results of the occultation analysis is reported in § [71 We 
discuss the broad band properties of the CXB and Earth spectra in § [HI Finally, the last 
section summarizes our findings. 
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Observations 



The Burst Aler t Telescope (BAT; iBarthelmy et al.l l2005l ). aboard the Swift mission 
(iGehrels et al.ll2004l ). launched by NASA in 2004, represents a major improvement in sen- 
sitivity for X-ray imaging of the hard X-ray sky. BAT is a coded mask telescope with a 
wide field of view (FOV, 120° x 90° partially coded) sensitive in the hard X-ray domain 
(15-200 keV). BAT's main purpose is to locate and to study Gamma-Ray Bursts (GRBs). 
While chasing new GRBs, BAT surveys the hard X-ray sky with an unprecedented sensitiv- 
ity. Thanks to its wide FOV and its pointing strategy, BAT monitors continuously a large 
fraction of the sky (up to 80%) every day. 

The Swift satellite constraints require that the pointing direction be at least 30° above 
the Earth's horizon. Nevertheless, due to its extent, it may happen that the Earth disk 
occults a substantial portion (up to 30%) of the BAT FOV. Moreover, BAT survey data 
include episodes of large occultation (up to ~70%) caused by the Earth when the spacecraft 
was in 'safe' mode. 

We use 8 m onths of BAT data which constitutes a well characterized dataset of BAT 
survey data (see lAiello et al.ll2008al jbl. for details) to study the different components of the 
BAT background. Our first aim is to derive the BAT background spectrum in the infinite- 
rigidity approximation. We then use all occultation episodes, as described in § [31 to derive 
a measurement of the CXB and the Earth atmosphere spectra. 



2.1. The BAT background 



The BAT background is highly complex and structured; it exhibits variability dependent 
on both orbital position and pointing direction. BAT employes a graded-Z fringe shield to 
suppress the in-orbit background. The fringe shield, located around and below the BAT 
detecto r plane, reduces the iso tropic cosmic diffuse flux and the anisotropic Earth albedo by 
~95% (IBarthelmy et al.ll2005l ). The two main background components are the CXB emission 
and the cosmic ray induced (prompt and delayed) backgrounds. 

The CXB spectrum in the 3-400 keV r ange is derived from H EAO-1 data. The following 
analytical approximation was suggested by lGruber et al.l (119991 ): 



Scxb{.E) 



7.877E-o-29e-^/4i-i3 

0.0259 (1)-^-^ + 0.504 (f)-^'^^ + 0.0288 (f)" 



3 < ^ < 60 keV 



1.05 



(1) 



E>m keV 

where Scxb{E) is expressed in units of keV cm"^ s^^ sr^^ keV^^. Given the large FOV 
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(~1.4sr half coded), the CXB is the dominant background component in BAT up to ~50- 
60keV. 

The prompt CR background is due to spallation effects of incident CRs on the material 
of the spacecraft; since the Earth magnetic field modulates the flux of incident CRs across 
the orbit, such background component is expected to vary with the cut-off rigidity Rc (i.e., 
the minimum momentum an incident charged particle must have in order to penetrate into 
the Earth's magnetosphere). The delayed component is caused by the excitation of the 
materials from the incident CR flux. This component builds up on times short compared to 
the relevant decay lifetimes, then varies as the slower of the irradiation or the lifetime. The 
Earth's magnetic field includes an indentation in the southern hemisphere called the South 
Atlantic Anomaly (SAA). During each SAA passage, BAT experiences a sharp increase in 
count rate due to the increase of the incident CR flux and a delayed background due to 
de-excitation of the spacecraft materials. 

In order to discriminate the various components of the BAT background, we correlated 
the BAT whole array rate (in each energy channel and normalized by the number of working 
detectors) with several orbital parameters. The final goal is to derive a "steady-state" BAT 
background model which is unaffected by orbital variations. 

2.1.1. Data Selections 

Our aim is to determine a rate-rigidity relation in order to extrapolate the BAT array 
rates to the infinite-rigidity case; this allows us to model the background variability due to 
the prompt and delayed CR components. 

First we selected the data excluding all observations where sources with signal-to-noise 
ratio (S/N) greater than 8 are detected. An 8cr source produces an increase in rate of less 
than 0.5% in a typical 300 s observatioii^; thus all point-like sources below this limit give a 
negligible contribution to the background level. 

The next step is to eliminate all observations whose exposure time is less than 300 s. As 
Fig. [T] (left panel) shows, for exposures below 300 s, the rates show a clear anti-correlation 
with exposure time, with an increase of a factor ~3 in the rates for few seconds of exposures. 
Exposures below 300 s are usually the result of a truncated observation because: 1) BAT 
detects a GRB or 2) BAT enters in the SAA and data acquisition is suspended. For any 



^BAT survey observation have typically an exposure of 300 s, although shorter and longer exposures might 
exist. 
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of these reasons, data of truncated exposures are excluded by the present analysis because 
they are not representative of the average BAT background. 




Fig. 1. — Left Panel: BAT rate at 100 keV versus exposure time. The "truncated" 
exposures below 300 s are noisy. Despite the impression most (>90%) of the observations 
have exposure larger than (or equal) 300 s. Right Panel: BAT rate at 100 keV versus time 
after each SAA passage. Note the sharp increase in rate and decay behaviour when the 
spacecraft exits the SAA. The second sharp peak at ~4000 s is due to a subsequent passage 
within the tail of the SAA. The rate modulation visible after ~6000 s is due to the Earth's 
magnetosphere. Negative times are used for exposure taken within the SAA. 

Whenever the spacecraft exits the SAA, BAT experiences a rate decline due to de- 
excitation of spacecraft materials. As shown in Fig. [1] (right panel) the rates reach their 
normal level after ~5600 s after each SAA passage. By excluding all observations taken 
within ~5600 s of an SAA passage, we thus eliminate short-lived radioactivity effects. 

The BAT rates also show a correlation with the angle between the Sun and the point- 
ing direction. This correlation becomes visible at angles >120° and decreases with energy, 
disappearing at 70 — 80 keV. The reason of this rate increase with the angle to the Sun is 
unclear. However, since the number of these observations is small (~5%), we decided to 
exclude them from the present analysis. 

The BAT effective area declines with energy and at ~ 200 keV reduces to 1/5 of its peak 
value at 50 keV; moreover, at these energies the fringe shield becomes partially transparent. 
Thus, it is possible to use the high energy channels as a "particle" detector to monitor the 
background level of the instrument (i.e. these channels do not yield much information on any 
celestial signal). We found that imposing that the rate of the last energy channel (194 keV- 
6.5 MeV) be in the range 10-20 ct s~^ eliminates roughly 1% of the observations which are 
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outliers in all the correlations we have studied. 



After these cuts, we find, as expected, that the rates in each energy channel decrease as a 
function of the cut-off rigidity Re- This effect is shown for two representative energy channels 
in Fig. [21 We model this behavior with an exponential and a constant {Rate = C + Be'^^^'), 
and fit this model to each energy channel. The fitted constant C provides an estimate of 
the BAT rate in the infinite-rigidity extrapolation. The distribution of the steepness of the 
rate increase with rigidity {a values), shown in left panel of Fig. El has a mean of -0 .34 and 
a RMS of 0.04, in perfect agreement with previous measurements ( llmhof et al.lll976l ). 
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Fig. 2. — Left Panel: Correlation of BAT rate with the cut-off rigidity at 100 keV. The 
dashed line is the best fit using an exponential plus a constant. Only those observations which 
passed the selection criteria explained in § 12.1.11 were used. Right Panel: Correlation of 
BAT rate with the cut-off rigidity at 170 keV. The dashed line is a exponential plus a constant 
fit. The outliers present in both figures around 12 GV are due to the spacecraft being in the 
vicinity of the SAA. Only those observations which passed the selection criteria explained in 
§ 12.1.11 were used. 

The BAT spectrum obtained by extrapolating the rates of each energy channel to infinite 
rigidity is shown in the right panel of Fig. [31 The bumpiness between 60 and 100 keV is due to 



the numerous fiuorescence emission lines from the fringe shield (see IWillisI [20021 . for details) 



3. Earth occultation 



The Swift orbital constraints require that the BAT pointing direction be always at least 
30° away from the Earth horizon. This is because the Earth is bright in Optical and X-rays, 
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Fig. 3. — Left Panel: D istribution of the a values. The mean of -0.34 is in good agreement 
with measurements from llmhof et al.l (Il976l ). Right Panel: BAT background spectrum 
extrapolated in the infinite rigidity case. 



thus it may damage the UV/Optical telescope (UVOT) and the X-ray telescope (XRT). On 
May 31, June 12 and July 28, 2005, the Swift spacecraft entered into 'safe'-mode because 
of star tracker loss of lock. In safe-mode operations, the XRT and UVOT telescopes are 
closed, but BAT still takes data. The spacecraft remains in sun reference pointing until 
commanding from the ground recovers Swift back to its normal status. In the timespan 
between the safe-mode and the recovering operation, the satellite uses the magnetometers 
and the sun sensor to derive its pointing direction. At least in the occasions mentioned 
abov^, the Earth passed through the BAT field of view (FOV). Fig. H] shows the BAT 
pointing directions during the deep occultation episodes described here. 

As the Earth partially occults the FOV, BAT registers a sharp decrease in rate due to 
occultation of the CXB emission. This is especially evident below 40keV where the CXB 
radiation dominates above the Earth's atmospheric components. The left panel of Fig. [5] 
clearly shows the drop in rates caused by the occulting Earth; between 18-20 keV the rates 
drop by a factor 3.5 when ~ 60% of the BAT FOV is occulted (as shown in the right panel 
of [5]). Thus, the Earth occultation can be used to measure the CXB emission by means 
of the depression caused in the BAT rates. Unfortunately, the Earth is not only a passive 
occulter, but also an active emitter. The Earth is a pow erful source o f X- and gamrna-ray s 
due to cosmic ray bombardment of its atmosphere (see iPetryl l2005l : ISazonov et al.l 120071 ). 



^ A few episodes of Earth occultation were found in BAT data, but some did not pass the criteria explained 
in Sec. [2XT] 
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Fig. 4. — Pointing directions for the occulted observations. Marker sizes vary linearly with 
the occulted FOV fraction. Occultations vary from 0.1% to ~80 %. In order to avoid 
contamination of the CXB signal from the Galactic Ridge emission, we used only occultation 
episodes at |b| >20° (data points already excised from this plot). 

This radiation is usually referred to as albedo, and it is discussed briefly in the next section. 



3.1. Atmospheric Albedo Gamma-rays 

The atmospheric Albedo flux is produced by cosmic ray interactions in the Earth's 
atmosphere. Hadronic interaction with atmospheric nuclei of the incident cosmic rays leads 
to the production of an electromagnetic and nuclear cascade with muons, nuclear fragments, 
and other hadrons. Gamma-rays above 50 MeV are produced mainly by the decay of mesons. 
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Fig. 5. — Left Panel: BAT rate, in the 18-20 keV range, as a function of the elevation angle 
above the Earth horizon. For ELV <30-60°, BAT starts to experience Earth occultations. 
The decrease in rate is expected to vary linearly with the occulted solid angle if Earth 
emission is negligible. Note that this graph is not a light curve as each data point is a 
separate (in time) 300 s observation. Right Panel: Example of deep Earth occultation of 
the BAT FOV. The black area is the region of the BAT FOV which is completely occulted 
by the Earth during the 300 s observation; the lighter gray is the un-occulted part of the 
FOV, while the region in between the black and the lighter gray is partially occulted due to 
the spacecraft movement in the 300 s. The colorbar shows the fractional time a given sky 
pixel is unocculted. 



while at X-ray energies the main source can be attributed to bremsstrahlung from secondary 
electrons. 



Measurements of the X-ray albedo radiation are repor ted in lSchwartz fc Peterson! (119741 ) 



in the 1-100 keV e nergy range and by 



Im 



hof et al.l (1l976l ) above 40keV. The Albedo spec- 



trum measured by ISchwartz fc Peterson! (11974! ) shows a cut-off below 30 keV, probably due 
to self-absorption of the radiation emitted from the inner layers of the atmosphere and a 
progressive flattening around 40 keV. Above 40 keV the Albedo emission decreases as a power 
law with p hoton index of ~ 1 . 4- 1.7 . This power-l aw behavior is confirmed by other experi- 
ments (e.g. ISchonfelder et al.!ll980l : lGehrelslll992! ) and by a recent Monte C arlo simulation 
of the hard X-ray emission of the Earth's atmosphere (jSazonov et al.! 12007! ). However the 
absolute normalization of the observed Earth spectrum depends on the altitude and the 
inclination of the satellite's orbit and on the solar cycle. 
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3.2. Method of analysis 

The rate R{E)i measured at energy E by BAT in a given observation, during which the 
Earth is in the FOV, can be described as: 

R{E)i = I{E)i — Vti ■ [R{E)cXB,i — R{E)Earth,i] (2) 

where the subscript i refers to the i — th observation, VLi is the "effective" solid angle 
occulted by the Earth, I{E) is the total background, and R{E)cxb and R{E) Earth are the 
CXB and the Earth emission respectively. 

The observations we are dealing with are generally non- contiguous, and thus all changes 
in the instrument configuration (e.g. number of working detectors) must be taken into 
account. We do this by computing the "effective" solid angle occulted by the Earth for each 
observation. This is defined as: 

n, = J2^y{i-ATf''n-v; (3) 

i=o 

where A''^ is the total number of sky pixels, AT^^^'^'^ is the fractional exposure tim^ a sky 
pixel of solid angle is unocculted and VJ is the vignetting affecting that sky pixel during 
the i-th observation. 

Equation [2] shows the "degeneracy" problem which limits the Earth occultation tech- 
nique when used to determine the CXB emission. Indeed, the measured depression of the 
rates with respect to the normal sky intensity level are a measurement of the difference of the 
CXB intensity and the Earth's atmospheric emission. Following the notation of Equation 



[21 this can be express ed as R{E)i = R 



E)cx B.i — R{E)Earth,i- adop t here an approach 



^t 

similar to the one of IChurazov et al.l (l2007l ) and iFrontera et al.l (120071 ). which consists of 



deriving the "difference" ( ON - OFF) spectrum and fitting it with a-priori spectral models. 

The difference spectrum is derived fitting Equation[2]to each energy channel (an example 
is shown in Fig.E]). In all these independent fits the two parameters {I{E)i and R{E)) are left 
unconstrained. Moreover, in order to avoid contamination by the Galactic Ridge emission we 
used only occultation episodes at Galactic latitude larger than 20°. The difference spectrum 
is shown in Fig. [71 However, before describing the spectral fit we discuss in detail the sources 
of systematic uncertainties affecting our analysis. 



■^The fractional exposure time is the fraction of the exposure time the sky pixel is unocculted. Thus, it 
varies from to 1 for completely occulted and unocculted pixels respectively. 
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Fig. 6. — Fit example for the 18-20 keV energy channel. The long dashed line is the best fit 
to the data. Note that the occultation of the CXB emission produces a very strong signal 
reducing the background rate by ~75%. 



3.2.1. A Note on the "Degeneracy" Problem 

The "degeneracy" problem (i.e. the fact that the CXB and the albedo emissions leave 
a similar signature during the occultation of the BAT FOV) might in principle be alleviated 
modeling the albedo emission of the Earth. This involves modeling the emission as a function 
of the cut-off rigidity of the visible disk as well as the (reasonably) expected non-uniformity 
of the albedo emission (i.e. limb or disk brightening effects). 

Indeed, since the albedo emission is generated at different cut-off rigidities, with respect 
the local rigidity of the satellite, one might reasonably expect that patches of the disk located 
at lower cut-off rigidities emit a larger X-ray flux. This information might be used to disen- 
tangle the albedo from the CXB signal. Moreover, the Earth is kno wn to be a non-uniform 



emitter at MeV and GeV energies. Both COMPTEL and EGRET fISchonfelder et al.lll980 
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Fig. 7. — The ON- OFF difference spectrum obtained fitting Equation [2] to each energy 
channel as shown in Fig. O The inset shows the negative part of the spectrum. Above 
100 keV the albedo spectrum dominates the CXB emission. 



Petryll2005l ) have shown that the Earth exhibits a bright limb. Thus, also this information 



might be used to model the expedected emission. 

However, we note that a few factors limit the application, in this analysis, of the model- 
ing described above. The limitations come from the fact that this analysis is entirely based 
on survey data. As explained in § I2.1.H the typical integration time for the survey is 300 s. 
Thus, we do not have a time-resolved monitoring (e.g. 1 s time resolution) of the transit of 
the Earth across the BAT FOV, but only 300 s snapshot observations with different level of 
occultations. Moreover, since the Earth is moving in the FOV within these 300 s, all the 
physical quantities (e.g. cut-off rigidity, fraction of the occulted FOV, etc.) are necessarily 
averaged over this time. Another limitation is due to the fact that the observations used 
here are not contiguous in time, but separated by weeks or months. Thus, the changing 
background conditions limit the precision of this analysis (as also shown in § 14. ip . 
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These facts limit the usage of a precise modehng of the Earth albedo emission which 
would allow to disentangle the albedo and the CXB signal without assuming a-priori spectral 
templates. A dedicated Earth observation with BAT in 'burst' mode (i.e. event-by-event 
mode) would not only allow to overcome the problems shown above, but also would extend 
the energy range of the measurement up to 350 keV (instead of 200 keV) and would also 
reduce the systematic uncertainties of the measurement to those related to the instrumental 
response only (see § H]) . 

4. Analysis of the Uncertainties 

4.1. Rate Variation 

The rate-rigidity graphs (examples are shown in Fig. [2]) show a scatter in the rate 
around the best fit which is generally larger than the statistical errors. This scatter is due 
to unknown effects. The pointing directions, the solar cycle, the spacecraft orientation with 
respect to the Earth and the Sun could be at the origin of this scatter which has an amplitude 
of less than 10%. We modeled the scatters as a Gaussian distribution such that the 1 o width 
of this distribution gives for each energy channel an estimate of the total (statistical plus 
systematic) error of the extrapolated rates. This constitutes the baseline uncertainty of this 
analysis, and it is propagated throughout all the further steps. 

4.2. Uncertainties connected to imprecise attitude determination 

During safe-mode operations, attitude determination relies on the magnetometers and 
Sun sensor. The derived attitude solution has a precision of the order of ~degree. This is 
confirmed by the analysis of sources detected during safe-mode pointings, which shows that 
the attitude differs from the nominal pointing direction by 1-2 degrees. The effective solid 
angle, computed in Eq. [31 is a slowly-varying function per degree of occultation. As shown 
in the right panel of Fig. [HI the fractional effective solid angle can be approximated by a 
straight line with a slope of 0.010 deg~^ in the 0.4-0.8 range of fractional occulted FOV. 
This means that an error of (at most) 2 degrees in the attitude determination translates 
into an uncertainty of ~2% in the determination of the occulted portion of BAT FOV. This 
additional systematic uncertainty is taken into account in our analysis. 
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Fig. 8. — Left Panel: Variation of the fractional effective solid angle as a function of 
elevation (distance of the pointing direction with respect the Earth limb). Right Panel: 
Fractional 1 a total uncertainty as a function of energy. The fractional uncertainty includes 
all the error estimates outlined in § |H The total error reaches its minimum of 4-5 % at 
the peak of the CXB emission. The peak around 100 keV is artificial and corresponds to 
the change of sign of the difference spectrum. Uncertainties above 100 keV are ~20 % and 
primarily systematic in origin. 



4.3. Uncertainties of the BAT instrumental response 

The BAT is a well calibrated instrument. However, given the very large FOV and 
the uncertainty in the modeling of spacecraft materials, the Crab Nebula is detected with 
slightly different spectral parameters across the FOV. To cope with this uncertainty, users are 
encouraged, when performing spectral fitting, to use a vector of energy-dependent systematic 
error^ which allows a unique spectral fit to the Crab Nebula wherever in the FOV. In this 
analysis, we account for such systematic errors, which fortunately have their minimum (~4%) 
in the 20-80 keV band. 



5. Instrumental response to a diffuse source 



The BAT response was developed by characterizing individual CdZnTe detector pixels, 
and by modeling the absorption and modulation of the coded aperture niask, then finally 
verifying by Monte Carlo simulations with radio active sources (jSato et al.ll2005l ). However, 



"'A detailed discussion is reported in |http: / / swift.gsfc.nasa.gov/docs /swift / analysis/bat-digest.html 
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since there remained uncertainty in response to continuum emissions, the response was ad- 
justed to fit the Crab nebula spectrum^. The BAT Crab spectrum can be described as: 
dN/dE = 10.40 E~'^'^'^ photons cm~^ s~^ keV~^. The values of normalization and pho- 
ton index are well within those used by most of the X-ray missions in a similar band (see 



Kirsch et al.ll2005l . for a review of Crab Nebula spectral parameters). 



However, the analysis of a diffuse source (as the CXB) presents some differences with 
respect to the study of point-like objects. Indeed, the spectrum of a point-like source is 
modulated by the coded mask pattern. Thus, the indirect unmodulated component which 
is scattered by the materials of the BAT instrument and of the satellite can be eliminated. 
Accordingly, the official response generator, batdrmgen, part of the standard BAT software, 
produces a response only for the direct component. However, the CXB, the subject of this 
paper, is seen as a diffuse emission and cannot be modulated by the coded mask pattern. 
We therefore utilized the Monte Carlo simulator to generate a more accurate response for a 
diffuse emission taking into account the scattered componenllfl as well as the effect of isotropic 
illumination of the BAT instrument. This simulator is the one used to verify the response 
on the ground, but the same corrections to fit the Crab spectrum are also applied. 



6. Results of the analysis 

In this section we present the main results of the analysis: the CXB and the Earth's 
atmosphere spectra. All the uncertainties described in § H] were added in quadrature to 
form the total uncertainty. The dependence of the total uncertainty with energy is shown 
in Figure [H (right panel). The total uncertainty reaches its minimum value of 4-5% at the 
peak of the CXB spectrum. 

In this section, all quoted errors on spectral parameters are 90% confidence for one 
interesting parameter. 



6.1. Spectral fitting 



The ON- OFF difference spectrum is folded in XSPEC (jArnaudlll996l ) with the proper 
instrumental response for a diffuse source. The model we used for the fit is the difference 
between the CXB and the albedo spectra. For the CXB spectrum we employ Equation [TJ 



^The so called 'un-coded' (or scattered) component comprises all those events which scatter in the satellite 
structure and produce a detectable signal in the BAT array. 



- 16 - 



For the albedo spectrum we use a jointly smoothed double power-law of the form: 



dN 



C 



dE {E/E,Y- + {E/E,Y^ 



[photons cm ^s ""^ sr ^] 



(4) 



where Fi and F2 are the two spectral indices and i?b is the break energy. This functional 
form reproduces well the atmospheric component with its dec linde at low energy, bump 
around 30-40 keV and a hard spectral index at high er energies (ISchwartz &: Peterson! Il974j : 



Imhof et al. 1976 



Earth emission (jSazonov et al. 



Gehrelsl 



1985 



Frontera et al.ll2007l ). Recent Monte Carlo simulations of the 



20071 ) show that Equation H] is a very good approximation 



of the Earth emission below 300 keV. We fix the values of the spectral indices and break 
energy at those suggested by ISazonov et al.l (120071 ). i.e. Fi=-5 and F2=1.4 and £^b=44keV. 
Thus, free parameters of our first fit are only the normalizations of Equations [1] and HJ 
respectively. The fit is poor, however, with a of ~ 220 for 75 degrees of freedom. Adding 
free parameters for the high energy spectral index F2 and the break energy E\) improves 
the fit (Xre(i=121.8/73). The F-test confirms that the improvement is very significant (F- 
test probability of 4.8xl0~^°). Adding another free parameter for the low energy spectral 
index Fi of the albedo does not improve the fit. Indeed, below 40keV the spectrum is 
completely dominated by the CXB emission, and thus it is not possible to constrain this 
parameter. Choosing as free variables the parameters of the E arth emission instead of those 



of the CXB spectrum is well motivated: there are indications (ISchwartzlll969l : lGehrelslll985 



Frontera et al.l 120071 ) that the high-energy spectral index of the albedo emission might be 
steeper than the cl assical value of 1.4. On the other side, the formula shown in Eq. [1] 
(iGruber et al.lll999l ) is a good representation of the broad band CXB spectrum. 



Our best-fit parameters (with 90 % CL errors) for the albedo spectrum are: F2=1.72±0.08, 
Eb=33.7±3.5keV and C=1.48+o:3 x 10"^. The normalization, and its 90% CL error, of the 
CXB as measured by BAT with respect Equation [1] is 1.06±0.08. This error includes also 
the change by ±1 in the low-energy spectral index (Fi) of the albedo emission. The CXB 
intensity in the 20-50 keV band is 6.43(±0.20)xl0~^ erg cm~^ s~^ sr~^. Figure shows the 
best fit and its residuals to the difference spectrum. 



7. Alternative measurement of the X-ray background spectrum 

The BAT in-flight background has a peculiar spatial distribution which shows larger 
count rates towards the center of the detector array and smaller rates towards its edges. 
This is clearly shown in Fig. (TU] (left panel). In the process of forming sky images, the 
BAT softwarqj removes this background component by means of an empirical bi-dimensional 
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Difference Spectrum 




20 50 100 



Energy (keV) 

Fig. 9. — Best fit to the ON- OFF difference spectrum. The model used is the difference 
between Equations [T] and HI Around 100 keV, the data become negative (as shown in 

Fig. ED. 

second-order polynomial function. 

This feature has an important physical meaning. Indeed, it is the result of a diffuse 
isotropic source (namely the CXB), shining through the mask. The peculiar shape of this 
"background" component arises from the fact that detectors at the edge of the array have 
a smaller solid angle of the sky as seen through the mask (they see the mask under large 
angles) than those at the very center. Given the extent of the BAT mask and array, this 
effect is noticeable and significant. 

We thus built a simple model, assigning to each detector its geometrical solid angle 



^For reference sec the description of the batclean tool available at 
http://heasarc.nasa.gov/lheasoft/ftools/headas/batclean.html, 
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Fig. 10. — Left Panel:BAT detector plane image in the 15-55 keV energy channel. Right 
PanehModel of the detectors solid angles distribution. The similarity of the model with the 
real data (left panel) is apparent. 



through the transparent mask elements. This model is shown in Fig. [TOl We can detect the 
CXB emission by fitting this model to the spatial distribution of the counts in each energy 
channel. However, this approach is valid only as long as the graded-Z shield remains opaque 
to X-ray photons(~50keV). Indeed, as soon as the shield becomes partially transparent, the 
effective detector solid angle increases because of the shield transparency. Thus, our model 
becomes inadequate above this energy. 



7.1. Model fitting 

Among all BAT observations which satisfied the selection criteria outlined in § 12.1.11 
we selected only those ones which were unocculted by the Earth. We then summed all the 
detector plane histogram^ (DPHs) into a single DPH with an overall exposure of ~1.8Ms. 
Summing the DPHs of observations with different pointing directions achieves the goal of 
smearing the contribution of sources which are below the detection threshold {8 a in this 
case). To each energy channel we fitted a model which is composed of: 



^BAT survey data are in the form of 80 channels detector plane histograms with a typical exposure time 
of 300 s. 
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• a constant term for the edges of the detector modules which register a higher count 
rate because of the larger exposed area, 

• the solid angle distribution model which takes into account the diffuse flux as seen 
through the mask, 

• a constant term which takes into account all other background components including 
the CR component which penetrates through the shielding. 



For each energy channel, the fit independently determines the intensity of the diffuse 
model. Since all the energy-dependent effects (e.g. absorption through the mask structure 
and transparency of the lead tiles) are correctly taken care of in the instrumental response 
described in § El we normalize our model (dividing by the maximum detector solid angle 
~0.6sr) and treat the dispersion of the solid angle distribution (~0.034sr) as a fluctuation. 
In this way, we make our approach insensitive to the exact computation of the solid-angle for 
each detector and at the same time, it allows us to use the same response matrix developed 
for the occultation measurement. We remark that this response matrix is based on extensive 
Monte Carlo simulations. 



7.2. Results 



For spectral fitting we convolved the CXB count rate spectrum with the BAT response 
matrix. For each energy channel, we summed in quadrature statistical uncertainty, the 
uncertainty on the mean solid angle (see § 17. ip and the uncertainty due to the BAT response 
(see §15]). 

A fit to the data (shown in Fig. [TTj) allowing only the overall CXB normalization (Equa- 
tion [T]) to vary yields a of 18.8 for 1 7 degree of freedom. The normalization with respect 
to the level of the CXB as measured by ICruber et al.l (Il999l ) is 1.091;q q3. This measurement 
is in very good agreement with the occultation measurement as Fig. [12] shows. 



8. Comparison with previous measurements 

In this section, we compare the CXB and the albedo spectra with previously available 
measurements in the same or overlapping energy bands. For reference, the values of CXB 
and albedo emissions as derived by BAT are reported in Table [T] 
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Table 1. Cosmic X-ray Background and albedo emission intensities 



Energy 


AE 


CXB 


CXB 1 a Error 


Albedo^ 


Albedo 1 a Errc 


keV 




keV^ cm~^s~-^ sr~-^ keV~^ 




10"^ ph cm-^s-i sr-i keV-^ 




15.0 


2.0 


43.92 


12.71 






17.0 


2.0 


46.82 


10.35 






19.0 


2.0 


48.58 


8.65 






21.0 


2.0 


50.39 


7.42 






24.0 


4.0 


50.96 


5.28 






28.1 


4.1 


47.12 


2.90 






32.2 


4.1 


45.26 


2.72 


63.96 


19.36 


36.2 


4.1 


45.07 


2.52 


73.93 


13.75 


41.4 


6.3 


44.00 


2.34 


76.25 


9.50 


47.8 


6.5 


40.39 


2.07 


74.94 


6.32 


54.3 


6.5 


37.33 


1.89 


64.57 


4.72 


62.0 


8.8 


33.25 


1.63 


55.01 


3.29 


70.9 


9.0 


32.99 


1.79 


38.30 


2.85 


81.2 


11.6 


28.08 


2.17 


34.16 


3.14 


92.9 


11.8 


26.25 


3.50 


26.53 


4.00 


106.1 


14.5 


24.86 


1.61 


20.64 


1.36 


120.8 


15.0 


23.18 


1.93 


16.41 


1.37 


137.2 


17.9 


21.15 


2.32 


13.43 


1.46 


156.8 


21.2 


18.95 


2.62 


10.96 


1.43 


181.1 


27.5 


18.02 


2.96 


8.38 


1.22 



^In the 14-30 keV energy range the \a upper hmit to the albedo intensity is 6.3x10 ^ ph cm ^ s ^ s 
keV-^ 
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Fig. 11. — Best fit to the second CXB measurement using Equation [T] with only the overall 
normalization left as free parameter. 



8.1. The X-ray Background Spectrum 

Both measurements of the CXB spectrum presented here produce the same results 
(within errors) for the normalization of the CXB intensity at its peak. Combining both 
measurements we determine that the CXB intensity at its peak is 8(±3) % larger than 



previously measured by HEAO-1 (IGruber et al.l 119991 ). We find that the CXB intensity in 



the 20-50 keV band is 6.50±0.15 x 10~^erg cm~^ s~^ sr^^. The observed intensity near the 
peak of the CXB spectrum (expressed in uF^, units) at 30keV is 46.2 keV^ cm~^ s~^ keV~^ 
sr~^. Figure US] shows the comparison of the BAT CXB spectrum with all other measurements 
available above 20 keV. All measurements agree well within 10 %. The detailed comparison 
is reported in Tab. [21 It is clear that the scatter in CXB intensities does not depend solely 
on the adopted spectra for the Crab Nebula. Some of the measurements showed in Tab. [2] 
might still be affected by systematic uncertainty in the instrumental response used. To our 
knowledge, BAT is the only instrument for which a dedicated instrumental response has 
been derived and tested for the analysis of the CXB. 
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Table 2. Comparison with previous results 



Instrument 


Ref. 


Energy 


band 


TP a 
-T Crab 


T b 
J-CXB 






keV 






HEAO-1/A2+A4 


1 


20- 


-50 


9.92^ 


6.06±0.06 


HEA0-1/A2 


2 


20- 


-50 


NA 


5.60±0.30 


BeppoSAX 


3 


20- 


-50 


9.22 


5.89±0.19 


INTEGRAL 


4 


20- 


50 


10.4 


~6.66 ^ 


BAT 


this work 


20- 


-50 


9.42^ 


6.50±0.15 



'^Crab flux quoted by the authors expressed in 10 ^ erg cm" 
keV-i 



■2 s-i 



^Intensity of the CXB quoted by the authors in 10 ^keV^ cm 



s ""^ sr ^ keV ^ 



iGruber et al.l (119991 ) do not report about their adopted Crab spec- 
trum; however the HEA0-A4 spectrum of the Crab Nebula can be 
described ( below 57 ke V) as dN/dE = 8.76 ^E-^.o.ts photons cm~^ 
s"^ keV"^ fl,Tunglll989l ). 

"^Authors do not give an exact measurement of the CXB flux, but 



report that their measurement is ~10% higher than the lGruber et al 
(119991 ) spectrum. 



'^The value quoted here has to be taken as a reference value. The 
systematic uncertainties discussed in § 14.31 allow to derive consistent 
Crab Nebula fluxes across the entire BAT FOV. 



References. — ( I jlGruber et all (119991): (2) 
(3) iFrontera et all (l2007h : f4) 



Churazov et al 



Marshall et al.l (1 19801 ): 



20071). 
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Fig. 12. — The two independent measurements of the CXB spectrum performed by BAT. 
The occultation measurement (gray datapoints) and the measurement derived using the sohd 
angle distribution (black datapoints) are in very good agreement. 



Fig. [13] shows a compilation of the X- and gamma-ray diff use backgrounds from keV to 

GeV energies. In addi tion to the work of iGruber et al.l (jl999l ). we show SM M (MeV) data 

(jWatanabe et al.lll997l). a nd COMPTEL and EGRET daia. in a recent revision (jWeidenspointner et al, 
IStrong et al.l 120041) In part i cular, the new analysis of EGRET data shows that the 



validity range of the iGruber et al.l (119991 ) formula (Eq. [T]) is now restricted to 3 keV< E < 
IMeV. 

We find that a good description of the available data in the 2 keV - 2 MeV range is 
achieved using a smoothly-joined double power-law of the form: 



E' 



dN 
'dE 



E' 



C 



[keV^photons cm ^s ^sr ^keV 



{E/EbY^ + {E/E^Y^ 



(5) 
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HEAO-1 -Gruberetal. 1999 
INTEGRAL - Churazov et al. 2007 
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BeppoSAX - Frontera et al. 2007 
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J I I I L_ 
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Energy [keV] 
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Fig. 13. — Comparison of CXB measurements above 20keV. The BAT spectrum (in red) is in 
agreement with the HEAO-1 (gray), INTEGRAL (blue) and BeppoSAX (black) observations. 
For clarity only the BAT occultation measurement is reported. 

The best fit, shown in Fig.[T4l yields values (and 1 a errors) of: C = (10.15±0.80) x 10"^, 
Fi = 1.32 ±0.018, F2 = 2.88 ±0.015 and Eb = 29.99 ± 1.1 keV. The reduced is acceptable 
(~1.2) considering the number (10) of different datasets fitted. The suggested formula 
reproduces well the CXB spectrum over two decades in flux and five in energy. At a given 
energy, the systematic uncertainty produced by the scatter of the measurements used here 
is of the of the order of 10 %. 

Note that, there is no astroph ysical need to conn ect the keV and the GeV diffuse back- 



grounds with a single formula (e.g. iGruber et al.lll999l ). It is generally agreed that the source 



populations contributing to the two diffuse backgrounds are probably different. Almost all 
of the CXB radiation up to 300 keV can be explained in terms of emission-line AGN (e.g. 
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Gilli et al.l 120071 ). Moreover, taking into acco unt (the like l y, but not yet detected) popula- 



tion of non-thermal electrons in AGN coronae, llnoue et al.l (120081 ) successfully reproduce the 
CXB emission up to 4 MeV. On the other hand, blazars account only for < 25 % of the GeV 
diffuse backgr ound and most likely other source classes contribute to the diffuse emission 
toermeil [20071 ). 
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1 — 
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+ 

o 

+ 

▲ 

+ 
+ 
+ 



I Mill I I I I 1 1 1 II I I I I 1 1 1 II I r 
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Nagoya balloon - Fukada et al. 1 975 — 
ASCA - Gendreau et al. 1995 
SMM - Watanabe et al. 1997 
GOMPTEL - Weldenspolntner et al. 2000 
EGRET - Strong et al. 2003 
INTEGRAL - Churazov et al. 2007 
RXTE - Revnivtsev et al. 2003 
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■ best fit to 2 keV < E < 2000 keV 
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Fig. 14. — BAT CXB spectrum compared with previous results. The dashed line is the best 
fit to 2keV< E < 2000 keV as reported in the text. 



8.2. The Earth albedo Spectrum 

Our Earth albedo spectrum is not compatible with the classical high-energy photon 
index of 1.4. The BAT data are consistent with a steeper high-energy photon index at 
99.989 confidence level. 
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Using the polar-orbiting satellite 1 972-076B. \lm]ioi et al.l (119761 ) found that above 40 keV 
the photon spectrum is consistent with a power-law with an index ranging from 1.34 to 1.4, 
depending on the latitude range scanned. Their measurement is based on the difference 
between pointings towards the atmosphere (Down) and pointing towards the sky [Up). In 
order to derive the alb edo spectrum, the authors sum the CXB emission and the Down- Up 
spectrum (see E g. 5 inllmhof et al.lll976l . for details). For the CXB emission, they adopt the 
measurement of lPall (119731 ). which describes the CXB photon spectrum as dN/dE = 25-E'^^'^. 
This representation differs from the HEAO-1 CXB spectrum in both normalizat ion and 
photon index in the 40-200 keV range. Thus, we adjust ed the llmhof et al.l (119761) albedo 
spectr a, taking into account the differences between the iGruber et al.l (119991 ) and the iPal 
(119731 ) CXB spectral representations. This is shown in Fig. [151 After the correction, the 
two (equatorial and polar) Albedo spectra are consistent with a power-law with photon 
index ~1.7. In particular the equatorial measurement is in good agreement with the BAT 
spectrum. 

The BeppoSAX satellite operated in a Low-Earth Orbit similar to Swift, but with dif- 
ferent inclination (4°). The BAT Earth albedo spectrum is compatible (w ithin the large 



uncert ainties of the BeppoSAX analysis) with the measurements obtained by iFrontera et al. 



(120071 ). It is worth noting that we derived the Earth intensity using different orbital positions 
(as done for BeppoSAX and OSO-3), and thus averaging over the magnetic latitude sampled 
by Swift. Moreover, as generally the Earth enters the FOV at large angles, we do not observe 
the upward albedo, but rather the albedo emerging at large zenith angles. This also seems 
confirmed by the similarity of our s pectruni with the downward gamma-ray flux measured 
for a balloon over Palestine, Texas (jGehrelslll985l ). 



Fig. [m reports also the predi ction of the Earth a lbedo emission as observed from the 



orbit of the INTEGRAL satellite (ISazonov et al. 



2007). Its normalization has been derived 



during the measurement of the CXB intensity (iChurazov et al.l 120071 ) . It is evident that 
this prediction and the BAT measurement agree well in shape but not so in normalization. 
Among many factors, the overall normalization depends strongly on the geomagnetic latitude 
and the distance to the Earth. The agreement of the BAT and INTEGRAL albedo spectra, 
respectively, with the equatorial and polar measurements of llmhof et al.l (119761 ) seem to 
confirm this interpretation. 



9. Discussion 



We have used Earth occultation episodes to derive with Swift/BAT an accurate mea- 
surement of the CXB emission in the 15-200 keV energy range. Moreover, we have proven 
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Fig. 15. — Measurements of the Earth albedo of llmhof et al.l (119761 ) at different geomagnetic 
latitudes. The original measurements (gray) have been corrected for un-subtracted CXB 
emission (see the text for details). After the correction, the albedo spectra (black data 
points) become steeper. 



by means of an independent technique the accuracy of the occultation analysis and of our 
results. The observed BAT intensity near the peak of the CXB spectrum at 30keV is 



46.2 keV^ cm ^ s ^ sr ^ keV ^ and its uncertaintjo is ~3% (including all systematics) . 
he normalization of the BAT CXB spectrum at 30 keV is ~ 8 % larger th an the HEAO-1 



Gruber et al.lllQQgi) measurement and consistent with the INTEGRAL one (jChurazov et al. 



20071 ). Moreo ver, considering that the precision of the HEAO-1 measurement at the CXB 
peak is 10% (jMarshall et al.lll980l ) and that BeppoSAX data are compatible with a larger 



Derived combining both measurements of the CXB. 
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Fig. 16. — BAT Earth spectrum as compared to past me a surern ents. References are shown 
in the legend. The datapoints (triangles) from iGehreld (119851 ) are a fit to t he downward 
gamm a-ray flux at 5g cm~^ over Palestine, Texas. The measurements from llmhof et al. 



(119761 ) were corrected to take into account the correct CXB emission (details are in the 
text). The thin solid line shows t he prediction of the Earth emission as obs erved from the 
orbit of the INTEGRAL satellite JSazonov et al.lboOTi : Ichurazov etaPboOTh . 



(up to 20 %) normalization of the CXB spectrum shows that all measurements above lOkeV 
are consistent within their systematic uncertainties. 



Such consistency is not observed at lower energies (e.g. see discussion in iRevnivtsev et al 



20051 ) . T he origin of this inc onsistency is unclear. However, it seems that neither cosmic 
variance jBarconset"aDl2000h nor differences in the flux scale calibration of each individ- 
ual instrument (IRevnivtsev et al.l l2005l : iFrontera et al.l 120071 ) may account for it. A likely 
reason for the discrepancy of CXB measurement in the 2-10 keV band might reside in a 
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systematic error in the response function used for diffuse sources (IFrontera et al.l 120071 ) . To 
our knowledge, BAT is the only instrument (beside HEAO-1 A2 which was designed with 
the purpose of measuring the CXB) which makes use of a dedicated instrumental response 
developed for this particular analysis. We also note that a rec ent measurement of the CXB 
performed, in the 2-7 keV, by Swift/XKT (jMoretti et al.l 120081) seems to confi rm the results 
of XMM-Newton dPe Luca fc Molen "dill2004h. RXTE jRevnivtsev et al.l boosh and Chandra 
( Hickox &: Markevitch]2006 ). If confirmed, this means that the CXB spectrum, as most re- 
cently measured, is 25-40 % larger (with respect the measurement of HEAO-1) below lOkeV 
while only ~10 % larger above 20keV. The functional form we provide in § 18.11 for the 
broad-band CXB emission approximates well this scenario. 

A larger, than previously estimated, CXB emission would in turn require a larger density 
of Compton-thick AGN both in the local and in the more distant Universe. Recently, Chandra 
stacking analyses of n iid-IR selected sources unveiled a la rge population of Compton-thick 
AGN at high redshift (jPaddi et al.ll2007l : iFiore et al.ll2008l ). This large fraction of Compton- 
thick AGN found at z=~ 2 can be accommodated if the emitted (obscured) fiux is very low. 
Indeed, lowering the assumed scattering efficiency (ratio of refiected to nuclear fiux) would 
increase by the same amount the number of Compton-thick AGN at any redshift. In this 
framework , the recent discoveries of Compton-thi ck AGN with an extremely low scattering 
efficiency (lUeda et al.l 120071 : IComastri et al.l 120071 ) fits well. These AGN are likely buried in 
a geometrically thick torus that obscures most of the nuclear fiux. Although there can be 
many of these hidden AGN, their individual contribution to the CXB is necessarily small. 
This seems to be confirmed by the fact that the contributio n of the mid-IR se lected, z=~ 2, 
AGN is < 3% of the CXB intensity in the 10-30 keV band JPaddi et al.ll2007[ l. Therefore, a 
larger contribution should be provided by Compton-thick AGN at lower redshift. 



10. Conclusions 



BAT performed a very sensitive measurement of the CXB emission in the 15-200 keV 
energy range. This measurement takes advantage of several episodes of CXB fiux modu- 
lation due to Earth's passages through the BAT FOV. We find that the BAT CXB spec- 
trum is in good agreement with the INTEGRAL one and that its normalization is ~ 8 % 
larger than the HEAO-1 measurement at 30keV. Additionally, performing an independent 
measurement of the CXB in the 15-50 keV band, we are able to confirm this result. Re- 
markably, our study also shows that all the available measurements in the > lOkeV range 
agree within their systematic uncertainties. The new analyses of COMPTEL and EGRET 
data (IWeidenspointner et al.l l2000l : IStrong et al.l 120041 ) show that the formula suggested by 
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Gruber et al.l (Il999l ) for the diffuse X- and gamma-ray backgrounds is only valid below 
2MeV. We derived a simple functional form which, in the 2-2000 keV range, approximate 
well (to a precision of 10 %) the CXB spectrum. 

Our study also derives the Earth albedo spectrum averaged over the magnetic latitudes 
sampled by Swift. The BAT spectrum is in agreement with all the previous observations 
performed by satellites operating in similar LEO orbits. This work shows that the Earth 
albedo spectrum declines at energies >40 keV according to a power-law with photon index of 
1.7, and n o t as 1 .4 as previously thought. A re-analyis of the measurements performed by 



Imhof et al.l (Il976l ) is in perfect agreement with the BAT Earth albedo spectrum. The good 
agreement among the available measurements allows to use the BAT Earth albedo spectrum 
to predict the background contribution from the Earth for other instruments operating at 
similar orbits. 
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